function [ tranProbMatrix ] = setInitialTransitionProbabilities(nstates)
    %Given the number of states, N, returns N*N matrix of transition
    %probabilities
    tranProbMatrix = zeros(nstates);
    for i = 1 : nstates-2
        tranProbMatrix(i,i) = 0.4;
        tranProbMatrix(i,i+1) = 0.4;
        tranProbMatrix(i,i+2) = 0.2;
    end
    
    tranProbMatrix(nstates-1,nstates-1) = 0.5;
    tranProbMatrix(nstates-1,nstates) = 0.5;
    
    tranProbMatrix(nstates,nstates) = 1;
    
    %Which yields this for N = 7
    %0.4 0.4 0.2 0 0 0 0
    %0 0.4 0.4 0.2 0 0 0
    %0 0 0.4 0.4 0.2 0 0
    %0 0 0 0.4 0.4 0.2 0
    %0 0 0 0 0.4 0.4 0.2
    %0 0 0 0 0  0.5  0.5
    %0 0 0 0 0  0    1.0
    
    
end

